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Abstract 

The path integral formulation of singular systems with second 
order Lagrangians is studied by using the canonical path integral 
formulation method. The path integral quantization of Podolsky 
electrodynamics is studied. 
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1 Introduction 



The study of constrained systems with higher order Lagrangians has been 
applied in many relevant physical problems. Poddolsky's electrodynamics 
[1] and a relativistic particle with curvature and torsion in three dimen- 
sional space-time [2] are some examples. 

The treatment for theories with higher order Lagrangians has been 
first developed by Ostrogradski [3 ]and leads to obtain the Euler and the 
Hamiltons equations of motion. 

The Lagrangian formulation of these theories require the configura- 
tion space formed by n generalized coordinates qi,q\ and q\. The Euler 
Lagrangian equations of motion, which are obtained from 

S = J L(qi,qi,qi)dt, (1) 

using the Hamiltons principle, are given by: 

dL d dL d 2 dL 

dq~ ~ JedqV + dP^Wi ' 

The passage from the Lagrangian approach to the Hamiltonian ap- 
proach is achieved by introducing the generalized momenta (pj, 7Tj) conju- 
gated to the generalized coordinates (g,, q\) respectively as 

dL d ,dL , . . 

% dt dqi 

dqi 

The phase-space will then be spanned by the canonical variables (qi,pi) 
and (<7i,7Tj), where = q\). 

However, a valid phase space is formed if the rank of the Hessian matrix 

OqWi »' J = 1 '-' n ' (5) 

is n. Systems which have this this property are called regular and their 
treatments are found in a standard mechanics books. Systems which have 
the rank less than n are called singular systems. 
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Now we will give two formulations to investigate singular systems which 
are: Dirac's method and the canonical path integral formulation [6-11] 

2 Dirac method 

The well-known method to investigate the Hamiltonian formulation of con- 
strained systems was initiated by Dirac [4,5]. In his formulation one defines 
the total Hamiltonian as 

H T = H c + v a H' a) a=l,...,m<2(n-l), (6) 

where H c being the canonical Hamiltonian and determined as 

H c = + ntf, - L, (7) 

v a are unknown coefficients. 

Due to the singular nature of the Hessian, we have a functionally in- 
dependent relations of the form 

H' a (qi,Pi,qi,ni) « 0. (8) 

Consistency conditions 

H> a = {H , a ,H e } + u lt {H' a ,H' fl }nO, (9) 

leads to the secondary constraints . Repeating this procedure as many 
times as needed, one arrives at a final set of constraints or / and spescifies 
some of v a . such constraints are divided into two types: first-class con- 
straints which have vanishing Poisson brackets with all other constraints 
and second-class constraints which have non-vanishing Poisson brackets. 
As there is an even number of second-class constraints, these can be used 
to eliminate conjugate pair of (p's, q's) and (it's, q's) from the theory by 
expressing them as functions of the remaining (p's, q's) and (it's, q's). The 
Dirac Hamiltonian for the remaining variables is then the canonical Hamil- 
tonian plus all the independent first class constraints So that the total 
Hamiltonian is defined as 

H T = H C + v x ^ x . (10) 
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Since first-class constraints are the generators of gauge transformations, 
this will lead to the gauge freedom. In other words, the equations of motion 
are still degenerate and depend on the functional arbitrariness. Besides, 
some v\ are still undermined. To remove this arbitrariness, one has to 
impose external gauge fixing conditions for each first class constraints. 

Fixing a gauge is not always an easy task, which make one be careful 
when applying Dirac's method. 

Now we would like to give the path integral formulation using the 
canonical path integral method and demonstrate the fact that the gauge 
fixing problem is solved naturally if this method is used. 

3 The canonical path integral formalism for 
second order Lagrangians 

Recently the canonical method [12-14] has been developed to investigate 
singular systems using the Caratheodoy's equivalent Lagrangian method 
and the equations of motion are obtained as total differential equations in 
many variables. 

Now we will give a brief review of the Caratheodory's equivalent La- 
grangian method. Let us consider a Lagrangian L(q iy q\, q\, t), we can ob- 
tain a completely equivalent one by 

L = L(qi,qi,qi,t) — , (11) 

such a function S(qi, qi, t) must satisfy 



H = Piqi + Triqi - L, (13) 
Pi = ^ (14) 
vr, = (15) 
These are the fundamental equations of equivalence Lagrangian method. 
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If the rank of the Hess matrix is n—R, R < n, then the generalized 
momenta conjugated to the generalized coordinates % are defined as 



dL . . 

n a = w^, a = R + l,...,n, (16) 
dq a 

dL . . 

n a = — a = l,...,R. (17) 
dq a 

Since the rank of the hess matrix Sn- is n — R, then one can solve 

dqidqj 

the n — R accelerations q a in terms of coordinates g,), the momenta n a 
and q a as follows 

c[a = w a {q u 7i a , q\). (18) 
Substituting (18) in (17) one has 

dL 



7T„ = 



dq a 



■- \qa=w a (q i ,q i ,iTa,q a ) ~ ^aili) Qii Vbi ^ a) ■ (19) 



On the other hand, from equation (3) if the rank of the Hess matrix 
gdjjl is n — r, we can obtain a similar expression for the momenta p 7 : 

p 1 = -HP(qi,q b ,p b ,7r a ), 7 = l,...,r, 6 = r + l,...,n. (20) 
The Hamiltonian if is defined as 



H = p b qb + ? 7 P 7 |p e =-tff + TTa^a (21) 

+q a -K a \ Pfj= _ H v-L(q h q h q\,q a = w a ), e = l,...,r. (22) 

Relabeling the coordinates t and g 7 as to and t 7 respectively, and 
will be called t a . defining the momenta Po as 

P. - f , (23) 

then the set of Hamilton Jacobi partial differential equations [HJPDE] is 
expressed as 
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dS dS 

Hq = P + H (t ,try,t a ;q b ,q a ;p b = — ; vr a = — ) = 0, (24) 

dS dS 

H ,p = p~, + H*(t , t 7 , t a ; q b , q a ] Pb = 7^-; n a = — ) = 0, (25) 

#a = ^a + HZ(to,t-y,t a ;qb,qa',Pb = ^~'^a = W~) = 0, (26) 

OQb CQa 

The equations of motion are obtained as total differential equations in 
many variables as follows: 



dH' dH' p dH' n 
dq t = -rr^dto + -rr^dt^ + -7r^dt a , (27) 
dpi dpi dpi 

dH' dH' p dH ,7T 

dq . = ^> dtQ + __L dt + —2.dt a , (28) 

07li U7li 07li 

f)H' dH' p f)H ,7T 

d Pl = ~^Pdt - —tdt, - ^dt a , (29) 
oqi dqi dqi 

dH' dH' p dH >7T 
dm = --^rdt - -^dt 7 - -^dta, (30) 
dqi dqi dqi 

f)ff' dH' p f)H ,7T 

dP ° = -^r dt ° - -sr 1 ^ ~ ^r" dt °" ( 31 ) 

OLq ULq uto 

1r7 . TT dH' dH' , 
dZ = (-H + p b — -1 + 7T a — ^ )dto 
opb dn a 

dH' p dH' p 
+ (-H p +p b —t + ir a —t)dt y 
7 dp b dn a 

dH' w dH' w 

The set of equations (27-32) is integrable if [14] 



dH'o = 0, (33) 
dH' p = 0, (34) 
dH'l = 0. (35) 
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If conditions (33-35) are not satisfied identically, one considers them as new- 
constraints and again tests the consistency conditions. Thus, repeating 
this procedure one may obtain a set of conditions. Hence, the canonical 
formulation leads us to obtain the set of canonical phase - space coordinates 
as follows 



Qb = Qb(to, tji ta), Pb = Pbito, ty, t a ), 

b — r + 1, n, 7 = 1, ...,r, a — 1, R. (36) 

Qa = Qa(to,t 7 ,t a ), 7T a = 7Y a (tQ,t^,t a ) , 

a — R + 1, n, 7 = 1, ...,r, a — 1, R. (37) 

Besides the canonical action integral is obtained in terms of the canoni- 
cal coordinates. H'q,H' p and H'^ can be interpreted as the infinitesimal 
generators of canonical transformations given by parameters t , £ 7 and t a 
respectively. In this case, the path integral representation may be written 
as [9-13] 



r ' R 

(q b , q a , try, t a \q' b) q' a , t' 7 , t' a ) = / ]J dq b dp b ]J dq a dir a x 

J 6=1 a=l 

exp i{ / (-H + Pbl -^ + n a —^)dt 

Jt y ,t a op b On a 



dH' p dH' p 

1 dp b dn a 
r)H' w r)H f7T 

+( - ff * + p '*f + T "*f )<i( ° } ' (38) 



4 Example 



In this section we will consider a singular system with a Lagrangian density 
depend on the dynamical field variables: C = C(if), d^ip, dfj,d u tp). One can 
obtain the Euler Lagrange equations of motion as follows 
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The momenta conjugated to ip and up are: 



p = T^r - 2d ^7^r0 ~ * = w, (40) 

dip o\OkW) oip 

7T — —^r. (41) 

Now will we will consider the Podolsky electrodynamics as a con- 
strained system with second order Lagrangian . The Lagrangian density 
for such a system is given as 

C = -\f^F, u + a 2 d x F a %F^ n, v,a,p = 0, 1, 2, 3, (42) 

where 

= d tx A v - d u A^ (43) 

and the metric convention g^ u = diag(+l, — 1, —1, —1). With the dynam- 
ical variables chosen as A 11 and A^ = the conjugated momenta are 
obtained as 



P„ = -Fo, - 2a\d k d x F 0X 5 l ; - d d x F*), (44) 
n, = 2a 2 (d x F 0X 5l-dxF^). (45) 

The primary constraints are 

H[ = n = 0, (46) 
H' 2 =p - d k 7r k = 0. (47) 

The expressible velocities A 1 are obtained as 

A^ = ^-^ + d k F' k + d l A . (48) 
2 or 

The canonical Hamiltonian is given by 

H = J d 3 x(p^ + ix^ - C) . (49) 



Making use of equation (48), we get 



H = J d 3 x[A d l ir t + p t A l + -^n i n%F ik + ^d k F lk 

+n^A + l -F, u F^ + \{A t - diAoXA* - <9M°) 

-a 2 (d k A* - d k d k A )(d t A* - dt&Ao)}. (50) 

Equations (46), (47) and (50) lead us to obtain the set of Hamilton 
Jacobi partial differential equations [HJPDE] as follows 



<99 

H'=P + H ; Po = ~^, (51) 

dS 

H[=n = 0, vro = — , (52) 

dS 

H' 2 =p -d k 7r k = 0. Po = ^s, (53) 

The equations of motion are obtained as total differential equations in 
many variables as follows 



dA l = 


dH' , 

dpi 


dHi , 7n dHL , <n 

dpi dpi 


(54) 




A'dt, 




(55) 


dA* = 


am , 

-TT A dt + 

OTTi 


T'dAV+^dA*, 

<77Tj dlTi 


(56) 






d k F* k + d*A ), 


(57) 


dp 1 = 


dH «dt 

dAi 


dA % dA dA ' 


(58) 




{-d l d k 7i k 


+ di&ir* - d k F kl )dt, 


(59) 


dp° = 


a a dt 
dA 


dA dA dA,' 


(60) 






- 2a 2 d i d i d k F k )dt, 


(61) 



9 



= (-p' - F oi - 2a 2 a i a t F t *)A, (63) 

= -^dt - §21^0 _ £S A4 °, (64) 
<x4 oA <x4 

= (-d k p k )dt, (65) 

rfp o = _M rft _ _ ?El dA o = Q (66) 

9* dt dt y J 

To check whether the set of equations (54-66) is integrable or not, let 
us consider the total variations of (51-53). In fact 

dH' = H' 3 dA = (d k Pk )dA = 0, (67) 
dH' = -H'gdt = {-d k p k )dt = 0, (68) 
dH' = H' 3 dt = (d k Pk )dt = 0. (69) 

The total variation of H 3 is identically zero. Hence, the equations of motion 
are integrable and the canonical phase space coordinates (A 1 ^ 1 , A 1 , 7r l ) are 
obtained in terms of parameters (t,A°,A°). Besides the canonical action 
integral is obtained in terms of the canonical variables as 

dz = J d^x^n^ - \f, u F^ - ^(A t - d t A )(A> - d l A°) 
-a 2 (d k A k - d k d k A )(d t A* - d t d l A )]dt. (70) 

Making use of equation (38) and equation (70), the path integral for 
the Podolsky electrodynamics is given as 



(A\ A\ t, A , A°\A i ', A\ t', A ', A ') = j ' f\dA l dp 1 dA* dn* x 

•* i=i 



-l(A i -d i A Q ){A*-#A Q ) 

-a 2 (d k A>*-d k d k A ) x 

(did - d t d l A )]dt}. (71) 
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The path integral representation (71) is an integration over the canonical 
phase-space coordinates (A\ p l , A\ 

5 conclusion 

We have obtained the path integral for singular systems with second order 
Lagrandians. For the Podolsky electrodynamics example since the integra- 
bility conditions oIHq = 0, dH[ = and alH' 2 = are satisfied the canonical 
phase-space coordinates (A l ,p l , A" 1 , it 1 ) are obtained in terms of parameters 
(t, A , A ) and the path integral (71) is obtained directly as an integration 
over the canonical phase-space coordinates without using any gauge fixing 
conditions. The Faddeev's Popov [15,16] method treatment for this model 
needs gauge fixing conditions to arrive at the result (71). The generaliza- 
tion of the present work for Lagrangians of order higher than two is given 
in reference [11]. 
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